Duality Relation for Quantum Ratchets 
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A duality relation between the long-time dynamics of a quantum Brownian particle in a tilted 
ratchet potential and a driven dissipative tight-binding model is reported. It relates a situation of 
weak dissipation in one model to strong dissipation in the other one, and vice versa. We apply 
this duality relation to investigate transport and rectification in ratchet potentials: From the linear 
mobility we infer ground-state delocahzation for weak dissipation. We report reversals induced by 
adiabatic driving and temperature in the ratchet current and its dependence on the potential shape. 
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Periodic structures with broken spatial symmetry, 
known as ratchet systems pj , present the attractive prop- 
erty of allowing transport under the influence of unbiased 
forces. The interplay of dissipative tunneling Q with un- 
biased driving enriches the quantum ratchet effect with 
features absent in its classical counterpart like, e.g., cur- 
rent reversals as a function of temperature [21 |4( . Quan- 
tum ratchet systems have only recently been experimen- 
tally realized in semiconductor Q and superconductor Q 
devices. Also from the theory side there are still few 

works HE0EHII3 which ' with the exception of 0,13, 
are restricted to the regime of moderate-to-strong damp- 
ing. After the pioneering semiclassical work 0, further 
progress towards a quantum description was made in , 
where the role of the band structure in ratchet poten- 
tials sustaining few bands below the barrier was investi- 
gated. Recently, a quantum Smoluchowski treatment |Hj 
added to the available methods. In this paper, we gener- 
alize to an arbitrary ratchet potential a duality relation 
put forward in 11] for a cosine potential. It provides a 
tight-binding description of quantum Brownian motion 
in a ratchet potential, and leads to an expression for the 
ratchet current valid in a wide parameter range including 
weak dissipation and nonlinear adiabatic driving. We ap- 
ply this method to discuss rectification and ground-state 
delocahzation occurring for weak dissipation in ratchet 
potentials. Our results encompass correctly the classical 
limit. 

We consider the Hamiltonian Hn of a quantum par- 
ticle of mass M in a one-dimensional periodic potential 
V(q + L) = V(q) tilted by a time-dependent force F(t), 
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Hn{t) ~ 2M 



V(q) - F(t)q. (1) 
The potential assumes in Fourier expansion the form 



EGO 
i= V l cos{2nlq/L-ip l ), 



(2) 



and can take any shape. Apart from special configu- 
rations {V;sin(iy9/ — lipi) = VZ} of the amplitudes Vi 
and phases tpi , this potential is spatially asymmetric and 
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FIG. 1: Dual relation between a dissipative ratchet system 
and a tight-binding (TB) model sketched for a two-harmonics 
ratchet potential (thick curve). Each harmonic (thin curves) 
generates couplings to different neighbors in the TB system, 
according to Eqs. JSJ and J7J|. The periodicity L of the TB 
model is determined by the viscosity n in the original model. 



describes a ratchet system. The interaction of the sys- 
tem with a dissipative thermal environment is modeled 
by the standard Hamiltonian Hq of a bath of harmonic 
oscillators whose coordinates are bilinearly coupled to 
the system coordinate q |2|. The bath is fully charac- 
terized by its spectral density J(oj). We consider strict 
Ohmic damping J(lo) = rjuj, which reduces to instan- 
taneous viscous damping (viscosity rf) in the classical 
limit. In such a system, the ratchet effect is characterized 
by a nonvanishing average stationary particle current 
= Hindoo t^ 1 J Q dt'v(t') in the presence of unbiased 

driving, characterized by lim.t->oo i -1 J Q dt'F{t') = 0, 
switched on at time t' = 0. In this paper, we shall 
consider the particular case of unbiased bistable driving 
switching adiabatically between the values ±F. We re- 
port a method to evaluate the stationary velocity v™ c (F) 
in the biased situation of time-independent driving F, 
which is also of experimental interest 0, 0| . The ratchet 
current in the presence of adiabatic bistable driving can 
be expressed as — v^ c (F) + v^ c (—F). 

The whole information on the system dynamics is con- 
tained in the reduced density matrix p(t) = TrelT(i), ob- 
tained from the density matrix W(t) of the system-plus- 
bath H = Hji + Hb, with time-independent driving F, 
by performing the trace over the bath degrees of freedom. 
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To evaluate the evolution of the average position (q(t)) — 
TrR,{q j 5(t)}, the diagonal elements P(q,t) = (q\p(t)\q) of 
the reduced density matrix are needed, and can be ob- 
tained by real-time path integrals techniques Q ■ The ve- 
locity follows by time differentiation. At initial time t' — 
0, we assume a preparation in a product form where the 
bath is in thermal equilibrium with the ratchet system 
W(0) = j o(0)e~' 3#B [TrBe-' 3AB ]- 1 . The bath tempera- 
ture is fixed by T = 1/ 0ks ■ This leads to a double path 
integral 

P(qt,t)= (3) 

9f qt 

J d Qi J d qil ( qi \p(0)h') JvqJ V*q'A[q}A*[q'}F[q,q'} 



on the continuous coordinates q and q'. Here A[q] = 
exp{—(it/h)HYi} is the propagator of the ratchet sys- 
tem for a path qit'), and F[q,q'] the Feynman- Vernon 
influence functional of the bath inducing nonlocal-in- 
time Gaussian correlations between the paths qit 1 ) and 
q'(t') 2\. Due to the nonlinearity of the potential V(q), 
these path integrals cannot be performed explicitly. For 
a cosine potential, Fisher and Zwerger introduced an 
exact expansion in the propagator A[q] which transforms 
the path integrals into Gaussian ones that can be per- 
formed. Generalizing this idea for the arbitrary periodic 
potential we find the expansion 
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where p{t') = {2-kTi/L) Y™ =1 ljS(t' - tj), and 

A ; = ^e^ ! for I > 0, A_, - A,*. (5) 

The physical meaning of these new quantities will be 
discussed later. For each term of the sum on m in Q 
we have introduced m intermediate times tj, and corre- 
sponding indices lj taking any value among {±1, ±2, . . .}. 
The sum ^2n.\ runs on all configurations of these in- 
dices. A similar expansion is performed for the propaga- 
tor A* [q'} , involving a new set of m' times i'- and indices l'j 
being used to define p'it') similarly to pit'). This enables 
us to rewrite the average position (qit)} — J dqqPiq,t) 
in terms of a series in the amplitudes Vi of the potential. 

Though still intricate, the resulting expression becomes 
easier to treat in the long-time limit we are interested 
in. Quantitatively, the measurement time t should be 
very long on the time scale 7 -1 = (r^/Af) -1 set by dis- 
sipation. A second approximation is necessary to pro- 
ceed: we neglect terms e _7tj , e _7t j , e~ 7 ( t-tj \ e _7 ^ _t j\ 
e - WB t 3j e -«B^ e -w B (t-t J ) ) an d e -WB(t-tJ) s w here lo b = 

/h, in the integrands involved in the series ex- 
pression for (qit)). We shall refer to this assumption 
as the rare transitions (RT) limit and discuss its valid- 
ity later. Generalizing we consider a Gaussian wave 
packet centered at position qo = Trn[qp(0)] and momen- 
tum po = Ttr[pp(0)] as initial preparation for the ratchet 
system. We obtain the important result 

(q(t)) ~ qo + - + — -(9tb(<))tb. (6) 

Parts of the series expression for (q(t)} has been summed, 
yielding the first three terms. The rest can be identified 



with the series expression for the expectation value of the 
position operator <7tb = L n \ n ) ( n l °f a driven 

tight-binding (TB) system, described by the Hamiltonian 

00 

Htb= £ (At\n + l){n\ + A.t\n){n+l\)-FfrB, (7) 

n,l — — 00 

and bilinearly coupled to a different bath of harmonic 
oscillators. The spectral density of this bath Jtb(^) = 
J(w)/[1 + (w/7) ] is still Ohmic but presents a cutoff at 
the frequency 7 set by dissipation. At initial time t' = 
the TB system is prepared in the state \n = 0). The 
calculation shows that the A; introduced in (JSJ are iden- 
tified with the couplings of the TB system Q . We stress 
that the Zth harmonic of the original potential results in 
a coupling to the Ith neighbors in the dual TB system as 
sketched in Fig. One can easily show that the spatial 
symmetry condition on the phases tpi is the same in both 
systems. The first three terms on the right-hand side 
of 10 reproduce exactly the classical solution for the av- 
erage position (q(t)) of a free system, V(q) = 0, at long 
times. In this linear case, the quantum and classical solu- 
tions should be identical, due to Ehrenfest theorem, and 
they are, because the TB average (<7tb(*))tb vanishes in 
the absence of the potential Viq). We expect the same re- 
sult when the potential is present but unimportant, e.g., 
for large driving F and/or high temperatures T. 

The series expression for the diagonal elements of the 
reduced density matrix /5tb of the TB system, which 
leads to the series expression for (<ztb(£))tb, can be 
written in terms of pairs of TB trajectories 9tb(^') = 
V' 1 Jo dt"p(t") [with pit') introduced above Eq. ©], 
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FIG. 2: Representation of some of the second-order (a,b) and 
third-order (c) paths contributing to the diagonal elements of 
the reduced density matrix of the tight-binding model, and 
the corresponding dependence on the couplings A;. 



and ?tb(^') defined similarly in terms of p'(t'). From 
that one extracts the spatial periodicity L of the TB sys- 
tem, yielding L — L/a, where a = TjL /2irh is the di- 
mensionless dissipation parameter of the original system. 
These pairs of trajectories combine in discrete paths in 
the q — q' plane parametrized by pairs of integers (n, n'). 
Each path starting in the diagonal element (0, 0) and 
ending at time t in (m,m) contributes to (m \pr&(t) \ m). 
Each transition in the path brings a corresponding fac- 
tor A; and all paths involve at least two transitions (cf. 
Fig. 0). Written in this form, the diagonal elements of 
the reduced density matrix are a solution of a general- 
ized master equation ^1 m terms of transition rates T m 
from the TB site (n, n) to the site (n + m,n + m). Con- 
sequently, these rates are expressed in power series of 
all the couplings A;, starting from second order. As the 
times tj , t'j introduced in Q are identified with the tran- 
sition times in the TB representation, the rates T m give 
also a way to control our assumption of rare transitions. 
It corresponds to neglect those paths which involve tran- 
sitions on a time scale max(7~ 1 , uo^ 1 ) after the initial 
time t 1 = or before the final time t' — t. As transitions 
in the TB model happen on a time scale T^ 1 , the du- 
ality relation will be valid when the transitions are rare 
on the time scale max(7 _1 , Wg 1 ), i.e., when all rates sat- 
isfy r m <C min(7,WB)- This condition is controlled by 
the dissipation through 7 = rj/M and the temperature 
through lub = litk^T /Ti. 

Due to the change of periodicity length between the 
two systems, the dissipation parameter a and the energy 
drop per unit cell e = FL become a = l/a and e = e/a 
in the TB system. Thus, weak dissipation in one system 
maps to strong dissipation in the other one although the 
viscosity 77 in the spectral density does not change. The 
asymptotic dynamics is usually described by the nonlin- 
ear mobility [i = lim^oo v(t)/F. With these notations, 
the duality relation JSJ) can be rewritten in the form 



(i(a, e) > (i - fiT B (l/a, e/a), 



RT 



(8) 



where /io = 1 / 77 is the mobility of the free system, V(q) = 
0. In the special case of a cosine potential, this relation 
was already obtained in for the dc mobility. It it also 



interesting to notice that it was also derived in [13j for 
the linear ac mobility in a cosine potential. However, we 
did not completely succeed in generalizing Eq. (jHJ in the 
presence of time-dependent driving. 

We shall now focus on the evaluation of the station- 
ary velocity v B ? c (F). By solving the generalized master 
equation mentioned above, one finds the stationary ve- 
locity v^ B — L J2 m m ^m in the dissipative TB system. 
The duality relation © can then be used to obtain 

mr m -r_ m . (9) 

m— 1 

As discussed above, the rates F m are power series in the 
couplings A; starting from second order. For a given m ^ 
0, there are only two possible second-order contributions 
to T m , which, after use of Eq. 10), sum up to 2] 
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d T e -m 2 aQ(T)+im(FL/hj)r 



(10) 



The influence of the dissipative environment enters 
through the dimensionless bath correlation function 
Q(r) = 2 J °° lj(1 ^ 2) [coth (w/20)(l - cos lot) + ismLor] 
with 9 = k^T/Ti^f. At zero bias F = and in the scal- 
ing limit 3> k^T, the rates show a power-law de- 
pendence on temperature r„ cx T 2m Q_1 . The linear 

(2) 

mobility /xtb is thus dominated by the rate T\ at low 
temperatures, and vanishes at T = for a < 1, which 
corresponds to free dynamics fi = /io in the dual weak- 
binding system Q. This suggests that the occurrence 
of a delocalization to localization transition at a = 1 for 
the ground state of a cosine potential 0, 0] would not 
be affected in more general potentials (see also Fig.|3. 

In the remainder of the paper, we focus on the 
ratchet current induced by adiabatic bistable driving 
v r = v dc(F) + v dc(~ F). The second-order rates obey 

(2) (2) 

r, n (— F) = r_^(F) an d therefore cancel out in the ex- 
pression for the ratchet current. Hence, we have to focus 
on contributions of at least third order to the rates T m . 
Here we neglect higher orders. This is known to provide a 
good approximation in TB systems with large dissipation 
parameter a and/or high temperature Q. For simplicity 
we also consider a potential consisting of only two har- 
monics. There is no problem of principle to include more 
harmonics 113. We find, with m = ±1, ±2, 



r (3) 
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where we have introduced cp = ip2 — 2<^i, and 
Jo 

e -2aQ(T+p)+aQ(T+2p) , g -2aQ(T- p)+&Q(r -2p) 



(11) 



Gf( 3 )r T ) = / dp e aQ(p)-2aQ(T+p/2)-2aQ(T-p/2)^ ^) 

Jo 



FIG. 3: Ratchet current and stationary velocity (inset) as a 
function of temperature for the potential of amplitude AV 
depicted in Fig. Q Weak dissipation is chosen with a = 0.2 
and fry = 0.76AV. Driving is set to FL = 0.57AV. 



At third order the rates obey r„ [F, ip) = r^(— F, —ip), 
which is a consequence of parity. The dependence of the 
ratchet current on the potential parameters is then up to 
third order in the potential amplitude 



oc ViV 2 sm((p 2 - 2^1). 



(13) 



The ratchet current vanishes for a symmetric potential 
sin (ip 2 — 2^1) = as it should. 

The behavior of the particle and ratchet currents as 
function of temperature and driving is shown in Figs. |3| 
andQJfor a two-harmonics potential. In Fig. EI the driv- 
ing is set to FL = 0.57Ay, whereas in Fig. 21 the tem- 
perature is fixed to k B T = 0.23A^. With V x = AV 2 , the 
untilted potential, depicted in Fig.^ has a barrier height 
AV = 2.2Vx. We choose a = 0.2 and fry = .76AT/. It 
means that the typical action is V2MAVL 2 w 2h, and 
the dissipation rate 7 = 77/M is about one-fourth of the 
classical oscillation frequency SIq = ^tt^Vi/ML 2 in the 
untilted potential (weak dissipation). In this numerical 
application, none of the rates exceeds 0.057 and 0.08ljb, 
which means that the duality relation is valid for this 
system. Moreover, the third-order rates stay at least one 
order of magnitude below the second-order ones. The 
ratchet current presents several reversals as a function 
both of the driving and the temperature. As expected 
for the small values of driving and dissipation used in 
Fig. |21 the stationary velocity is very close to the value 
of a free system vo — F/r) at T — 0, which corresponds 
to localization v^ B w in the TB system jj^- Accord- 
ingly, sa in this regime. The stationary velocity also 
tends to vq (dashed line in Fig. for driving or temper- 
atures much higher than the potential barrier, and the 
ratchet current vanishes correspondingly. If observed in 
experiments, this linear behavior would provide a direct 
estimation of dissipation. 

In conclusion, we obtained a duality relation yielding a 
tight-binding description of Brownian motion in a tilted 



FIG. 4: Stationary velocity and ratchet current (inset) as a 
function of driving for the potential of amplitude AV depicted 
in Fig. Q The dashed line is the classical solution in the 
absence of potential. Weak dissipation is chosen with a = 0.2 
and fry = 0.76AV. Temperature is set to k B T = 0.076AV. 



ratchet potential. We demonstrated its application to 
investigate rectification of adiabatic driving and ground- 
state derealization for weak dissipation. 
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was supported by the Dutch Foundation FOM. 
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